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HIGHLIGHTS 


•  Buckling  analysis  combined  with  experiments  enabled  to  fabricate  stable  YSZ  membranes. 

•  Effective  residual  stress  in  a  free-standing  YSZ  film  was  estimated  by  a  new  method. 

•  Comparisons  between  simulations  and  experiments  show  excellent  agreement. 

•  Novel  post-buckling  design  space  for  thin  electrolyte  fabrication  has  been  set  up. 
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The  buckling  behavior  of  a  thin-film  electrolyte  of  a  micro-solid  oxide  fuel  cell  is  investigated  based  on 
experimental  measurements,  analytical  estimations  and  numerical  simulations.  An  energy  minimization 
procedure  is  applied  in  combination  with  the  Rayleigh-Ritz  method  to  represent  the  buckling  modes, 
evaluate  the  buckling  amplitude  and  determine  the  threshold  values  for  instability  transitions  in  the 
system.  The  residual  stresses  in  the  film  deposited  on  a  silicon  substrate  are  evaluated  based  on  wafer 
curvature  whereby  nanoindentations  tests  are  applied  to  estimate  the  Young’s  modulus  of  the  deposited 
film.  The  partial  release  of  residual  stresses  in  the  film  during  free  etching  of  the  substrate  is  estimated  by 
a  new  method  combining  pre-etching  optical  measurements  with  posteriori  stress  analysis.  The  energy 
interpretation  of  the  obtained  deflection  shape  is  discussed.  Comparisons  between  simulation  results 
and  experimental  data  show  the  efficiency  of  this  method  to  predict  various  buckling  stages  of  free¬ 
standing  thin  films.  A  post-buckling  design  space  for  thin-film  electrolyte  fabrication  is  presented  by 
applying  a  stress-based  failure  criterion. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Micro-solid  oxide  fuel  cells  (/iSOFC)  are  currently  receiving 
increasing  attention  as  promising  electrochemical  energy  converter 
devices.  The  high  conversion  efficiency,  the  high  power  density,  the 
fast  start-up  times  and  the  low  operational  temperature  (350- 
550  °C)  make  this  a  suitable  power  source  for  small  portable 
electronic  applications.  Before  reaching  the  promised  potential  of 
this  technology,  some  challenges  need  to  be  overcome  like,  in 
particular,  the  thermomechanical  stability  during  fabrication  stages 
and  under  operational  conditions.  A  ^SOFC  is  built  of  a  multi¬ 
layered  structure  of  ceramic  and  metallic  plates.  Although  these 
structures  exhibit  advantageous  properties  such  as  a  favorable 
strength-to-weight  ratio,  a  high  fatigue  endurance,  and  a  strong 
corrosion  resistance,  they  may  exhibit  some  functional  problems. 
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For  example,  such  a  structure  has  a  low  resistance  towards 
peeling  stresses  at  the  interface  between  different  layers  where 
delamination  may  occur.  In  addition,  these  films  are  subjected  to 
in-plane  stresses  associated  with  both  the  thermal  expansion 
mismatch  and  the  fabrication  process.  This  may  induce  either 
membrane  buckling  [  1  ]  under  compression  or  membrane  fracture 
under  tension.  Both  buckling  and  cracks  reduce  the  load  carrying 
capacity  of  the  structure  [2].  Further,  crack  propagation  can  lead  to 
complete  failure  limiting  the  operating  time  of  the  fuel  cell. 

The  present  work  deals  with  the  mechanical  stress  analysis  and 
the  design  implications  for  the  fabrication  of  thin  films  in  ^SOFCs. 
The  fabrication  of  /iSOFC  membranes  involves  several  steps  of  film 
deposition  and  layer  etching  as  depicted  in  Fig.  1  see  Refs.  [3]  and 
[4].  After  the  etching  of  the  silicon  substrate  (Fig.  1  (c)),  the  mem¬ 
branes  are  free-standing  and  the  mechanical  system  in  Fig.  1(g)  is 
found  to  be  a  model  of  three-layer  clamped  plate  which  consists  of: 
a  top  electrode,  the  electrolyte  made  of  YSZ  (as  described  in  Ref.  [3]) 
and  a  bottom  electrode.  The  stability  of  such  a  system  should  be 
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Fig- 1-  The  fabrication  steps  of  /zSOFC  membranes,  (a)  Silicon  wafer  double  side  coated 
with  Si3N4.  (b)  Backside  photolithography,  and  reactive  ion  etching  of  Si3N4  (LPCVD). 
(c)  Wet  etching  of  Si  in  20%  aqueous  KOH.  (d)  YSZ  deposition  by  PLD  onto  free-standing 
Si3N4  membrane,  (e)  Reactive  ion  etching  of  Si3N4  membrane,  (f)  LSC  Cathode  depo¬ 
sition.  (g)  Anode  deposition  by  platinum  sputter-coating. 


guaranteed  under  both  fabrication  and  operation  conditions  at 
temperatures  between  400  and  500°  C. 


Fig.  2.  Schematic  representation  of  PLD  techniques  used  in  the  fabrication  of  thin  YSZ 
film  deposited  onto  a  silicon  substrate. 


the  fabrication  costs  of  the  thick  deposited  films  and  to  the  low 
electrochemical  performance  attributed  to  the  long  transport  path 
in  the  cross-plane  direction  of  thick  membranes. 

A  “design  accepting  buckling”  is  highly  important  for  the 
fabrication  of  thin  films  in  many  applications  where  a  stress-based 
failure  is  adopted.  Although  the  important  contribution  from  Ref. 
[9]  toward  a  design  in  post-buckling  regime,  the  application  of  the 
energy  method  in  such  a  work  is  restricted  to  the  axis-symmetry 
describing  the  first  buckling  where  the  secondary  buckling  (sec¬ 
ond  bifurcation)  is  not  captured  and  then  excluded  from  the  design 
space. 

Important  simulation  results  on  the  buckling  of  membrane  in 
^ctSOFC  are  shown  in  Refs.  [10],  however,  there  is  always  a  need  to 
provide  a  guideline  for  membrane  designs  based  on  both  simula¬ 
tion  results  and  experiments  where  fabrication  and  operation 
constraints  in  fuel  cells  are  considered. 

In  the  presented  study,  an  advanced  exploitation  of  the  energy 
method  is  contributed  to  provide  the  needed  informations  on  thin- 
film  design  in  different  buckling  stages.  The  numerical  results  are 


ll  Design  consideration  for  thin-film  manufacturing 

When  PLD  is  employed,  the  deposited  film  exhibits  different  re¬ 
sponses  with  respect  to  the  deposition  conditions  and  techniques.  A 
high-temperature  deposition  may  induce  a  residual  compressive 
stress  which  drives  the  membranes  to  buckle.  On  the  other  hand,  the 
deposition  at  room  temperature  may  initiate,  during  a  subsequent 
annealing  up  to  400°  C,  cracks  which  propagate  in  the  brittle  film 
causing  rupture  under  tensile  stresses  as  shown  in  Fig.  3  (see  also 
[3]).  The  residual  stress  is  then  a  finger  print  of  the  deposition’s 
mechanism  type  and  the  related  conditions  of  temperature  and 
pressure,  see  for  example  [5].  Another  examples  of  the  intrinsic  ef¬ 
fects  are  the  phase  transformation  inside  the  film,  the  densification, 
the  crystallization,  and,  the  chemical  association-dissociation  which 
induces  a  strain,  and  subsequently,  a  chemomechanical  stress  in  the 
material  ([6  and  [7]). 

In  this  work,  focus  is  made  on  the  stress  analysis  and  design  of 
YSZ  electrolyte  fabricated  by  pulsed  laser  deposition  PLD  (see 

Fig.  1(e)  and  Fig.  2). 

In  the  framework  of  ^SOFC,  some  design  concepts  are  adopting  a 
failure-based  buckling  criterion,  like  in  Refs.  [8],  where  the  use  of 
thick  film  is  recommended  to  avoid  both  buckling  and  fracture. 
Such  a  design  option  suffers  from  drawbacks  which  are  related  to 


Tension  Compression 


(a)  (b) 

Fig.  3.  Two  typical  results  of  free-standing  YSZ  membrane  obtained  from  pulsed  laser 
deposition,  (a):  YSZ  membrane  rupture  after  YSZ  deposition  at  room  temperature  and 
annealling  up  to  400  °C  at  3  °C  min-1,  (b):  membrane  buckling  after  PLD  deposition  at 
700  °C.  Images  are  taken  at  room  temperature. 
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experimentally  validated  in  order  to  present  a  usable  design  space 
for  YSZ  electrolyte  fabrication  when  the  material  properties  are 
given  and  the  residual  stresses  are  evaluated  depending  on  the 
manufacturing  conditions. 

By  considering  both  rupture  and  buckling,  the  design  concept 
adopted  in  this  work  allows  a  safe  and  flexible  choice  between 
various  options  of  deposition  conditions  and  film  dimensions. 
Knowing  that  multi-layer  thin  plate  clamped  onto  a  stiff  substrate  is 
a  prototype  model  for  several  laminated  devices,  the  design 
method  applied  here  is  exploitable  in  different  microfabricated 
products  using  thin  film  materials  such  as  in  chemical  sensors, 
piezoelectric  techniques,  purification  devices. 

2.  Energy  minimization  method 

When  the  total  potential  energy  of  the  thin  film  increases  under 
residual  compressive  stresses,  the  film  may  take  a  new  equilibrium 
state  by  performing  a  buckling  instability  transition  to  a  new  stable 
configuration.  Energy  minimization  method  is  an  adequate  nu¬ 
merical  approach  to  represent  the  principle  of  minimum  potential 
energy  and  to  predict  system  stability.  In  this  paper,  the  needed 
topics  on  the  simulations  method  are  shortly  introduced. 


2.1.  The  potential  energy  of  the  film 


The  elastic  energy  in  the  domain  Q  of  a  thin  film  with  clamped 
boundaries  (Fig.  4)  is  given  by 

F  =  IJ  ff :  affl.  (1) 

Q 

Here  a  is  the  stress  tensor  and  e  denotes  the  strain  tensor 
written  as 


Eij  (z)  —  <5ij  E0  +  4J1  +  Ejj  (z)  5 


(2) 


whereby  z,  j  e  {x,  y }  stand  for  the  in-plane  compounds,  Fig.  4.  The 
scalar  eo  is  the  average  value  of  the  residual  strain,  is  the  Kronecker 
delta.  The  terms  rjT  are  the  compounds  of  the  membrane  strain.  They 
describe  the  elongation/contraction  deformation  and  the  stretching 
by  shearing  deformation  of  the  middle  surface  of  the  film  when  any 
bending  resistance  is  excluded.  A  nonlinear  formulation  of  mem¬ 
brane  strain  is  used  in  this  based  on  the  von  Karman  plate  model: 


lyaui  9Uj  awaw\ 

«  2  \dj  9  i  dj  a  i  J  ’ 


(3) 


whereas  iZi  and  zi,  stand  for  the  in-plane  displacements  and  w 
denotes  the  out-of-plane  deflection  i.e.  u  =  (ux,  izy,  w). 

The  terms  in  Eq.  (2)  represent  the  bending  and  the  twisting 
deformation  of  the  film  and  is  given  as  linear  function  of  the  “z” 
coordinate  (see  Fig.  4): 


Fig.  4.  Schematic  representation  of  the  buckling  of  a  thin  clamped  film  with  side 
length  a.  The  vector  (ux,  uy,  w)  denotes  the  displacement  of  a  middle  plane  point. 


a2w 


(4) 


Considering  the  case  of  plane  stress  model,  the  elastic  energy  F 
introduced  in  Eq.  (1 ),  can  be  subdivided  into  the  membrane  energy 
Fm  and  the  bending  energy  Ft>: 

F  =  Fm+Fb.  (5) 


Using  the  in-plane  Hook’s  relations,  we  have  the  total  mem¬ 
brane  energy  Fm  of  the  film  domain  Q: 


a/2  a/2 

Fm  =  2(1  - v1)'  j  J  (exx  +  eyy  +  2i,£xx£yy)z=0dxd-)/ 

-a/2  -a/2 
a/2  a/2 

+  ff-v  J  J  £xy|z=0dxd3/-  (6) 

-a/2  -a/2 

The  terms  a  and  h  denote  the  film  side  length  and  thickness 
respectively. 

When  clamped  boundary  conditions  are  considered,  the  sum  of 
bending  energy  terms  affected  by  the  shear  modulus,  i.e.  F/2(l  +  v) 
(the  second  Lame  constant)  vanish  upon  integration  (this  includes, 
however,  the  “twist”  energy).  Hence,  the  remaining  bending  energy 
of  the  film  is  written  as: 


Fh  m 


Eh3 

24(1  —  v 2) 


a/2  a/2 

/  / 

-a/2  -a/2 


/  02W  02W 

^ax2  +  ay2 


dxdy. 


(7) 


The  numerical  solving  of  the  buckling  problem  consists  of  two 
steps:  firstly,  Rayleigh-Ritz  parameterization  is  applied  in  order  to 
represent  the  displacement  fields  as  summations  of  trial  functions 
weighted  by  unknown  coefficients  (amplitudes).  Secondly,  a  nu¬ 
merical  minimization  procedure  is  applied  in  order  to  determine 
the  unknown  coefficients  to  describe  the  film  configuration  in  the 
stable  equilibrium  state. 


2.2.  Rayleigh-Ritz  parameterization 


The  parameterization  is  split  into  two  summation  modes:  a 
mirror-symmetric  and  an  anti-mirror  symmetric  modes  with 
respect  to  oxz  and  oyz  planes  with  terms  denoted  by  the  super¬ 
scripts  “s”  and  “a”  respectively.  The  in-plane  displacement  along  x- 
direction  is  formulated  as: 


ux(x,  y) 


Y)  "ijAi**-  y)  +  Y  “fj/y  (y.  *)> 

i  j  =  1  i  j  =  1 


(8) 


where /-j  and/y(x,  y)  =  ff-(y,  x)  are  the  trial  functions  affected  by 
the  unknown  coefficients  uf-  and  uf-,  respectively: 

rs,  x  •  2z'7TX  (2j-l)7Ty 

Jij(x,y)  =  sin~^~  cos—  a  ■ 

The  in-plane  displacement  in  y-direction  is  written  as 


(9) 


uy(x,y)  =  ux(y,  -x).  (10) 

The  out-of-plane  deflection  is  parameterized  as  superposition  of 
buckle  functions  gf-  and  gf-  affected  by  the  unknown  amplitudes  wf- 
and  wf-: 

m  m 

w(x,  y)  =  wi#j(*>  y)  +  5/  wfj sfj(x,  y),  (11) 

i  =  l,j  =  i  i  =  lj  =  i+l 
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where 


s?j(x,  y) 


and 


2nrx  , 

cos - (-1)1 


cos^-(-iy 


cos 


2 'my 


COS 


a 

2 my 
a 


j^{0,e0)  denoted  by  ACl  =  Amin  =  rninjAi}™2!  vanishes.  Within  a 
range  of  high  in-plane  compressive  values  where  r0  <  $  <  there 
is  a  reduced  residual  strain  £q2  <  ,  for  which  the  Hessian  matrix 

at  the  equilibrium  point  (£Q2,T\/e  =  "0)  has  a  multiple  eigen¬ 
value  AC2  =  0.  The  strain  beyond  this  critical  value  (r0  <  £q)  aU°ws 
the  clamped  film  to  buckle  either  into  the  path  of  primary  buckling 
with  symmetric  mode  or  into  the  path  of  secondary  buckling  with 
symmetry-breaking  mode  denoted  BP2  in  Fig.  5. 


Sh(x,  y)  =  sin 


,.  (2l  +  1  )tcx 


/  iM  • 

-  (  -  1)  sin — 


sin(2j+V)ny 


-  (-l)jsin  — 


.  (2j+l)ra  ,  .  7rx\  (  .  (2i  +  l)7rv  ,  .  tuv\ 

-  (  sin v  J  J - (  -  l)Jsin —  •  sin- —  ;  -  (-l^sin  —  . 


It  is  observed  that  in-plane  and  out-of-plane  trial  functions 
satisfy  the  fixed  boundary  conditions  of  clamped  film. 

2.3.  First  and  second  buckling 

Using  the  vector  tf  of  the  in-plane  coefficients 
„  U  „  and  the  vector  ~W  of  out-of-plane  amplitudes 

.  U{w?}m  .  .  ,  the  potential  energy  in  Eq.  (1)  is  rein- 

L  u J  1=1  ,j=i  L  uji=ij=i+i  r 

troduced  here  as  a  differentiable  function  of  the  vector  V,  the 
vector  W ,  and,  the  scalar  qje  K  of  residual  strain: 

F  :  U2P2  x  IRm2  x  R - ^  IR 

(u,W,£o)^f(u,W,£o). 

At  equilibrium  state,  the  minimum  potential  energy  implies  an 
in-plane  stationary  condition: 


A  double  buckling  point  has  been  reported  in  several  cases,  see 
e.g.  Refs.  [11],  and  is  associated  with  an  eigenvalue  multiplicity.  In 
Ref.  [12]  a  method  for  the  prediction  of  second  buckling  is  introduced 
in  combination  with  application  of  the  finite  element  method. 

However,  within  the  presented  energy  approach,  the  buckling 
transition  of  the  square  film  can  be  iteratively  controlled  with 
respect  to  either  dimensional  or  loading  condition  conditions.  As 
example,  let  us  consider  the  case  of  a  film  of  Poisson  ratio  of 
v  =  0.20  (value  as  that  in  e.g.  Refs.  [13],  [14]  and  15])  and  thickness 
h  =  300  nm  subjected  to  a  residual  compressive  strain  qj  =  -0.002. 
The  simulations  were  performed  by  using  Ritz  expansions  of  degree 
m  =  4  and  p  =  5.  For  a  long  side  length  a,  this  film  can  be  in  a  stable 
secondary  buckling  mode  represented  by  the  path  PB2  in  Fig.  5 

where  the  eigenvalues  of  the  Hessian  matrix  34?(We,  qj),  Eq.  (13), 
are  strictly  positive  i.e.  Amin  >  0.  Such  a  secondary  buckling  mode 
was  iteratively  predicted  to  take  place  for  a>aCl  =  97  pm,  actually, 
the  det  e,  qj)  vanishes  by  an  iterative  decreasing  the  value  of 
the  film  side  length  a  toward  aC2  on  the  path  PB2.  On  the  other 


Z  -  0.  (12) 

dU 

The  in-plane  displacement  is  then  obtained  from  Eq.  (12)  as 
nonlinear  function  of  the  out-of-plane  deflection,  this  is  denoted  by 
the  relation  V  =  G(W ).  _ 

The  values  of  the  coefficients  wf-  and  wf-  (compounds  of  W )  are 
obtainable  by  applying  a  minimization  procedure  (conjugate 
gradient)  on  the  potential  energy  function  F(W,e0)  = 
F(G(W),  W,q,). 

For  a  given  residual  strain  q),  the  stability  of  the  system  is 
controlled  according  to  the  Hessian  matrix  providing  the  double 
derivative  of  the  energy  function  with  respect  to  deflection  am¬ 
plitudes  wfj  and  wf-.  At  the  equilibrium  point  (  W  e,  qj),  the  Hessian 
matrix  of  order  m2  is  obtained  as 

=  02F(wU) 

aw/2  ~w=~w( 

Let  us  consider  the  fundamental  path  of  trivial  solution 

W e  =  o'  which  is  denoted  as  FP1  in  Fig.  5  and  represents  the  state 
of  unbuckled  film.  The  onset  of  the  first  unstable  equilibrium 
(primary  buckling)  is  attributed  to  a  threshold  value  of  reduced 
compressive  strain  r0  =  Eq  <  0  when  the  minimum  eigenvalue  of 


Jt(We,E0) 


Fig.  5.  Schematic  of  bifurcation  of  the  equilibrium  solution:  (W e,E0).  FP1:  funda¬ 
mental  path  of  a  stable  unbuckled  state  (r0  >  fo  >wy  =  wfj  =  0).  FP2:  fundamental 
path  of  an  unstable  unbuckled  state  (f0  <#’  wfj  =  wfj  =  0).  BP1:  primary  bifurca¬ 
tion  path  with  axis  symmetry  shape  (f0  <  wfj^O,  wfj  =  0).  BP2:  secondary 
bifurcation  path  with  symmetry  breaking  shape  (r0  <  FqL  wfj  ^0,  wfj  =£0). 
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hand,  the  first  buckling  was  investigated  using  the  condition 

det  =  0,  e)  on  the  path  FP1  in  Fig.  5,  it  was  predicted  to  be 
taking  place  at  a  >  aCl  =12.8  pm. 

3.  Film  in  different  buckling  stages 

The  results  shown  in  Fig.  6  are  obtained  by  performing  simu¬ 
lation  tests  with  high  degree  of  Ritz  expansion:  m  =  10  and 
p  =  20  where  the  buckling  load  for  v  =  0.20  and  E  =  200  GPa  is 
g$  =  -0.538  MPa. 

The  first  buckling  mode  is  shown  in  Fig.  6(a)  for  a  residual 
compressive  stress  go  =  -12.5  MPa.  A  film  in  this  primary  buckling 
stage  allows  a  minimization  of  the  total  stored  potential  energy  by 
an  exchange  from  membrane  energy  to  bending  energy. 

Note  that  the  secondary  buckling  load  is  predicted  as 
(7q2  =  -30.9  MPa.  The  deformations  in  Fig.  6  (b  and  c)  represent 
an  advanced  stage  of  secondary  buckling  under  compressive  re¬ 
sidual  stresses:  go  =  -125  and  -325  MPa  respectively.  The  film’s 
shape  in  this  secondary  mode  is  characterized  by  the  symmetry 
breaking  associated  with  a  new  mode  of  energy  storage:  a  certain 
amount  of  membrane  energy  is  released  by  elongation  deforma¬ 
tion,  another  amount  is  transferred  to  bending  energy  and  a 
considerable  amount  of  membrane  energy  is  stored  as  a  membrane 
shear  energy.  Hence,  to  absorb  more  energy  in  bending  form,  the 
thin  film  should  have  more  local  curvature  and  this  justifies  the 
observed  branched  folds  and  wrinkling  at  the  boundaries.  The 
folds  in  Fig.  6(c)  are  more  twisted  compared  to  Fig.  6(b)  and  this 
associated  with  a  higher  amount  of  stored  membrane  shear  energy. 
See  Appendix  A. 

The  above  mentioned  Young’s  modulus  was  obtained  based  on 
nanoindentation  measurements  performed  on  an  8YSZ  film  of 
thickness  of  ft  =  300  nm  deposited  by  pulsed  laser  deposition  at 
700  °C.The  influences  of  the  variation  of  the  value  of  the  Young’s 
modulus  on  the  buckling  analysis  can  be  discussed  in  the  sense  of  both 
deflection  magnitude,  and  instability  transition.  Considering  a  film 
with  a  defined  ratio  of  thickness-to-side-length  (h/a),  a  given  Poisson’s 
ratio  {v\  and  exposed  to  a  constant  residual  compressive  stress  (oo): 

i  Within  a  given  buckling  mode,  a  changing  of  the  Young’s 
modulus  values  should  affect  the  magnitude  of  the  out-of-plane 
deflection:  an  increased/decreased  value  of  Young’s  modulus 
should  effectively  result  in  a  lower/higher  buckling  amplitude. 

ii  If  the  film  is  exposed  to  the  residual  compressive  values  (do) 
such  that,  the  mechanical  equilibrium  of  the  system  is  found  to 
be  in  the  neighboring  of  the  buckling  transition  (around  the 
point  of  an  unstable  equilibrium  state),  a  perturbation  of  the 
value  of  the  Young’s  modulus  may  shift  the  threshold  value  of 
the  buckling  loads  ( g q1  or  g\ f )  and  then  the  film  may  be  found 
within  another  mechanical  equilibrium,  i.e.  it  may  perform  a 
transition  between  unbuckled  state,  primary  or  secondary 
buckling  mode. 

On  the  other  hand,  if  the  film  was  exposed  to  a  constant  re¬ 
sidual  compressive  strain  (not  stress)  the  variation  of  the 
Young's  modulus  should  affect  only  the  buckling  amplitude 
without  changing  the  buckling  mode. 

4.  Experimental  validation 

We  shall  recall  that  for  a  relevant  comparison  of  simulation 
results  to  experimental  measurements  it  is  important  to  insure  the 
accuracy  of  the  used  input  data. 

For  this  purpose,  a  consistent  procedure  is  introduced  here  in 
order  to  extract  the  effective  value  of  residual  stresses  in  the  YSZ 
film  from  the  experiments. 


Mathematica  8  © 


(b) 


Mathematica  8  © 
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Fig.  6.  (a)  First  buckling  with  axis  symmetry  for  a0  =  -12.5  MPa.  (b)  and  (c)  Second 
buckling  with  rotational  symmetry  for  a0  =  -125  and  -325  MPa  respectively.  Ritz 
degrees:  m  =  10  and  p  =  20. 


4A.  A  sequence  of  curvature  measurements  of  multi-layers  sample 
without  etching 

The  residual  stress  in  the  multi-layers  domain  is  estimated  from 
the  knowledge  of  the  wafer  curvature  by  applying  of  Stoney  for¬ 
mula  (see  Ref.  [16]).  The  curvature  is  evaluated  at  room  tempera¬ 
ture  based  on  the  optical  measurements  sequentially  applied  on 
the  sample  at  different  coating  steps  depicted  in  Fig.  7: 

Firstly,  the  curvature  measurements  were  applied  on  a  top  side 
coated  silicon  substrate  and  then  on  a  double  side  coated  substrate. 
From  the  difference  between  (a)  and  (b)  curvature  measurements,  a 
tensile  stress  ^NS  =  128  MPa  is  estimated  in  the  silicon  nitride 
membrane.  The  stress  in  the  YSZ  layer  was  determined  by  taking 
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Fig.  7.  Optical  measurement  of  curvature  is  performed  at  each  step  of  sample  prepa¬ 
ration:  (a)  Silicon  wafer  380  pm  thick  is  top  side  coated  with  200  nm  Si3N4.  (b)  Wafer  is 
double  side  coated  with  200  nm  Si3N4.  (c)  Sample  is  top  side  coated  with  300  nm  YSZ 
film. 


Here  h\  =  300  nm  is  the  YSZ  film  thickness  and  hs  =  2h2  +  ^3  is 
the  substrate  thickness  with  h2  =  200  nm  and  h3  =  380  pm  are 
the  Si3N4  film  coating  thickness  and  the  wafer  Si  thickness 
respectively  (see  Fig.  8),  the  terms  B[  =  EfJ  -  vXy  i  =  1,  2,  are  the 
biaxial  modulus  of  YSZ  and  Si3N4  film  respectively.  The  Young’s 
modulus  of  the  YSZ  film  E  =  200  GPa  is  obtained  by  using 
nanoindentation  measurements  where  the  Poisson  ratio  is 
vi  =  0.20.  The  total  thickness  of  the  double  side  coated  substrate 
hs  =  2h2  +  /13  =  380.4  pm.  The  residual  stress  obtained  from 
curvature  measurements  in  substrate-supported  film  feature  are 
aj  =  -1350  MPa  for  YSZ  film  and  o^NS  =  128  MPa  in  the  Si3N4 
film.  The  biaxial  modulus  of  the  silicon  nitride  is  considered  as 
B2  =  385  GPa.  By  substituting  in  Eq.  (14),  the  effective  value  of 
compressive  residual  stress  in  the  free-standing  YSZ  film  is 
<tysz  =  —186  MPa. 


difference  between  (b)  and  (c)  curvature  measurements  when  this 
film  is  found  to  be  storing  a  residual  compressive  stress  of 
of  =  -1350  MPa. 

4.2.  Effective  residual  stress  in  free-standing  YSZ  film 


4.3.  Effect  of  underlying  SisN4  layer  on  membrane  buckling 

In  order  to  incorporate  the  effects  of  the  underlying  thin  Si3N4 
layer  (h2  =  200  nm)  on  the  buckling  of  YSZ  membrane,  Fig.  8,  a 
validated  procedure  is  introduced  based  on: 


It  is  worth  to  note  that  residual  stress  obtained  from  the 
aforementioned  measurements  does  not  represent  the  effective 
value  of  the  stresses  when  the  film  is  deposited  on  an  etched 
substrate.  The  etching  of  an  underlying  silicon  substrate  allows 
more  compliance  in  the  structure  since  the  total  thickness  has  been 
extremely  reduced  by  the  factor  (hi  +  h2)/(h^  +  2h2  +  *13)  <  1  / 500. 
Therefore,  the  stiffness  of  the  free-standing  (post-etching)  double¬ 
layer  is  extremely  lower  than  the  total  stiffness  of  the  substrate- 
film  sample,  Fig.  8.  This  reduces  an  important  amount  of  the  re¬ 
sidual  stress  especially  that  deformational  effects  arise  early  in  the 
structure  (even  before  buckling)  during  the  film’s  deposition  on  the 
thin  (etched)  substrate. 

Further,  it  is  not  practical  to  measure  the  pre-buckling  residual 
compressive  stress  in  a  film  deposited  on  an  etched  substrate  since 
we  take  out  the  film  already  buckled  from  the  deposition  machine. 

For  an  experimental  validation,  it  is  instructive  to  establish  a 
method  for  an  accurate  estimation  of  the  residual  stress  in  the  free¬ 
standing  membrane  deposited  on  an  etched  substrate. 

In  such  a  perspective,  the  stress  analysis  of  the  curvature  mea¬ 
surements  in  Fig.  7  is  applied  to  predict  the  misfit  strain  between 
YSZ  and  Si3N4  layers  in  substrate-supported  film  case.  Therefore, 
within  the  feature  of  free-standing  membrane,  a  stress  analysis  of 
the  YSZ-Si3N4  double  layer  is  reperformed  by  using  the  afore¬ 
mentioned  misfit  strain  as  interface  condition  and  finally,  the 
compressive  stress  in  the  free-standing  YSZ  film  is  obtained  as  an 
explicit  function  of  the  measured  stresses  and  cr|NS,  see  Appendix 
B  and  Appendix  C,  mainly  Eq.  (C.7): 


"ySZ 


ft  ~4,hs)af  ~  aSNS 

h\  ,  B2  \  3/i1(/ii+h2)2B2 
h2  '  Bi  '  h\Bx+h\B2 


(14) 


a.  The  driving  force  for  the  buckling  is  the  residual  compressive 
stress  restricted  to  YSZ  membrane. 

b.  An  underlying  thin  ShN4  layer  may  affect  the  plate  deformation 
by  increasing  the  stiffness  of  the  system. 

Therefore,  an  equivalent  buckling  problem  is  considered  as 
follows: 

i  A  plate  domain  is  restricted  to  the  clamped  YSZ  membrane  of 
thickness  h\  subjected  to  a  compressive  residual  stress  (7Ysz- 

ii  Increased  stiffness  values  equivalent  to  those  of  YSZ-Si3N4 
double  layers  are  assigned  to  the  single  YSZ  layer. 

Indeed,  the  prefactors  of  the  integrals  in  the  equations  (5),  (6) 
and  (7)  are  D  =  £/i3/12(1  -  v 2),  Cs  «  Eh/ 1  +  v  and  Cm  =  Eh /l  -  v2 
respectively.  They  represent  the  bending  stiffness,  the  membrane 
shear  stiffness  and  the  elongation  stiffness  respectively,  for  a  single 
layer.  They  are  substituted  by  assigned  values  of  average  total 
stiffness  for  double-layer  membrane: 


•  D 


Ai2(1-^)E2\ 
\h,(  \-v\)E, ) 


•  Gs 


Gs(l  + 


h2(l+Vi  )E2\ 


•  Gm 


Gm 


fr2(l-^)£2\ 


Note  that  the  subscripts  i  =  1,  2  stand  for  YSZ  and  Si3N4 
respectively.  The  above  assigned  values  are  obtainable  considering 
the  shear  modulus  and  plane  strain  modulus  as  average  values  over 


the  total  thickness  i.e.  Ya= 


1 2(ft1+h2)(l+i'i) 


and  Ya= 


1  (hi  +h2)(l-*f  )* 


hi 

h2 
hs 
h2 

Substrate-supported  YSZ  film 
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Free-standing  YSZ-Si3N4  membrane 


Fig.  8.  Residual  stresses  are  strongly  released  in  free-standing  sample  due  to  “etching-induced  compliance”  associated  with  an  extreme  decrease  of  total  thickness  by  a  factor 
(hi  +  h2)l(hi  +  2h2  +  h3)  =  1/761.4. 
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Simulations  of  the  buckling  of  YSZ  film  were  performed  with  both 
cases  of  single-layer  values  and  assigned  double-layers  stiffness 
values.  The  effect  of  the  underlying  Si3N4  on  the  buckling  results  was 
found  in  agreement  to  the  prediction  contributed  elsewhere  using 
another  approach  of  averaging  residual  strain  in  bi-layer,  see  Refs. 
[12],  page  111.  Moreover,  it  is  stated  experimentally,  that  reactive 
ionic  etching  of  the  thin  underlying  Si3N4  (see  Fig.  1  (e))  was  applied 
without  observing  considerable  change  in  the  buckled  film  config¬ 
uration.  Consequently,  the  deformation  of  the  double-layer  was 
mainly  driven  (dominated)  by  the  buckling  of  the  YSZ  membrane. 

4.4.  Optical  profilometry  measurements 

The  300  nm  thick  YSZ  membrane  was  deposited  by  using  PLD  at 
700  °C  onto  a  free-standing  layer  of  390  pm  side  length.  This 
membrane  is  found  to  have  a  buckled  shape  and  deformation  is 
measured  at  room  temperature  by  using  an  optical  profilometry 
technique  (Wyko  NT100  white  light  interferometer).  The  buckling 
shape  and  amplitude  is  detected  with  high  accuracy  and  illustrated 
in  3D  view  in  Fig.  9(a)  where  the  secondary  buckling  mode  is  clearly 
observed.  The  measured  amplitude  at  the  center  of  the  buckled 
membrane  is  Wmax  =  5.74  pm. 

On  the  other  hand,  numerical  simulation  is  performed  for  a 
residual  compression  00  based  on  the  aforementioned  analysis  of 
curvature  measurements  i.e.  cr0  =  0ysz  =  -186  MPa.  By  choosing  a 
degree  of  Rayleigh-Ritz  expansion  with  m  =  10  and  p  =  20,  the  use 
of  in-plane  equilibrium  conditions  reduce  the  degree  of  freedom  of 


(a) 


Fig.  9.  (a)  3D  view  of  a  buckled  300  nm  thick  YSZ  membrane  taken  with  white  light 
interferometry  shows  an  “octopus”-like  pattern  with  folded  boundaries,  (b)  Simulation 
results  showing  the  deflection  values  along  the  diagonal  direction  of  the  buckled  film 
(in-house  Mathematica  code  is  used). 


Table  1 

Buckling  amplitudes  obtained  from  simulations  are  compared  with  the  optical 
measurements.  Deviations  (AW)  are  given  relatively  to  the  maximum  measured 
deflection  Wmax  =  5.74  pm. 


Model  [pm] 

Measured  [pm] 

AW 

Upward 

5.47 

5.74 

4.7% 

Downward 

-0.51 

-0.8 

5% 

the  problem  to  m2  =  100.  The  simulation  results  in  Fig.  9(b)  pre¬ 
sents  a  buckling  amplitude  in  good  agreement  with  the  measured 
amplitude  value  of  buckled  pattern  shown  in  Fig.  9(a)  (see  Table  1 ). 

The  profile  of  the  measured  sample  deviates  slightly  from  the 
assumed  rotational-symmetry  and  this  reflects  a  local  imperfection 
in  the  buckling  behavior  of  the  pattern  related  to  uncontrollable 
material  and  geometrical  effects  during  sample  fabrication. 

5.  A  guideline  for  design  space 

For  a  reliable  structural  design,  it  is  important  to  define  the 
failure  mode  and  to  establish  the  safety  assessment  with  respect  to 
the  suitable  yielding  criterion.  The  maximum  stress  criterion 
(principal  stress  or  Rankine  criterion)  states  that  failure  occurs  at  a 
point  in  brittle  material  when  a  principal  stress  (03  or  031)  reaches 
either  the  uniaxial  tension  strength  at  >  0,  or  the  uniaxial 
compression  strength  ac  <  0. 

Flowever,  a  more  conservative  prediction,  Coulomb-Mohr  cri¬ 
terion,  is  often  used  to  predict  failure  in  brittle  materials  especially 
for  plane  stress  problem: 

(max{ff1,(7[1})+  |  (miniffi.ffn})- >  1 
fft  ffc 

Flere,  (.)+/(.)_  is  the  positive/negative  part  notation  (it  vanishes 
if  (.)  is  a  negative/positive  value,  else,  it  is  equal  to  (.) ).  Fig.  10  shows 
three  buckling  test  cases  simulated  with  a  low  order  Ritz  expansion 
m  =  4,  p  =  5.  Under  compressive  stress  value  a0  =  -300  MPa,  it  is 
observed  that  thin  film  with  a  short  side  length  (70  pm)  exhibits  a 
primary  buckling  mode  where  the  film  of  wide  side  length  (150  pm 
and  500  pm)  are  in  secondary  buckling. 

The  film  in  primary  buckling  undergoes  a  high  tensile  stress  in 
Fig.  10(bi)  near  the  boundaries  where  yielding  are  detected  in 
Fig.  10(di)  for  mechanical  strengths  at  =  250  MPa  (value  as  that  in 
e.g.  Refs.  [23]  and  [9])  and  ac  =  -1500  MPa.  On  the  other  hand, 
yielding  is  restricted  to  limited  zones  of  high  tensile  stress  in  sec¬ 
ondary  buckling,  see  Fig.  10(b2),  (d2).  Moreover  with  a  film  of 
500  pm  side  length  we  observe  more  stress  relaxation  within  an 
advanced  stage  of  secondary  buckling  where  the  film  is  fully  safe. 

The  stress  relaxation  in  post-buckling  stage  is  investigated  in 
Fig.  11  for  the  same  aforementioned  Ritz  degrees  and  material 
properties.  The  extremum  stresses  are  plotted  against  different 
values  of  side  length  a.  It  is  observed  that  stresses  are  relaxed 
beyond  a  safety  threshold  value  of  a  =  180  pm. 

This  simulation  has  been  reproduced  with  different  values  of 
residual  stresses  and  a  least  square  method  has  been  applied  to 
formulate  the  dependency  between  residual  stresses  and  safety 
lengths.  Subsequently,  a  design  space  is  formulated  with  respect  to 
the  geometrical  and  the  material  properties  of  the  film  and  shown 
in  Fig.  12.  The  yielding  regions  above  the  tensile  strength  and  bel¬ 
low  compressive  strength  are  excluded  from  safety-design  area.  In 
the  remaining  part,  three  safety  regions  are  located:  S4  which  is  the 
unbuckled  region  and  is  located  above  the  primary  buckling  curve 
which  is  a  region  of  low  stress  and  represents  the  classical  con¬ 
servative  area  of  budding-based  failure  design.  The  other  regions, 
Si  and  S2,  are  located  above  the  curves  Ci  and  c2  in  post-buckling 
regions  where  the  stresses  are  controlled  under  the  tensile 
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Fig.  10.  Simulations  are  performed  with  Ritz  degrees  m  =  4  and  p  =  5  on  a  film  of  thickness  h  =  300  nm  for  a  residual  stress  a0  =  -300  MPa.  The  subscripts  i  e  {1,  2,  3}  stand  for 
films  with  side  length  70  pm,  150  pm  and  500  pm.  (a,)  film  buckling.  Results  are  given  on  the  lower  surface  denoted  as  (b,)  and  (q)  for  principal  stresses  ai  and  au  and  (d,-)  for  safety- 
yielding  parts.  E  =  240  GPa,  v  =  0.20  and  <7t  =  250  MPa. 


strength  at  =  250  MPa.  The  area  of  safe  regions  Si  and  S2  are  clearly 
larger  than  S4,  therefore,  the  post-buckling  analysis  has  allowed  the 
prediction  of  a  new  wide  design  area. 

It  is  worth  to  note  that  our  experimental  tests  show  a  high 
survival  rate  of  square  8YSZ  films,  (samples  of  300  nm  thickness 


Fig.  11.  Extremum  values  of  principal  stresses  for  a  300  nm  thick  film  in  first  and 
second  buckling  stages  under  a  residual  stress  a0  =  -300  MPa.  Safety  lengths 
a  >  180  pm  are  determined  using  Rankine  criterion. 


with  different  side-length  values),  in  post-buckling  regime  when 
PLD  deposition  was  performed  at  different  temperatures  700  °C 
under  20  mTorr.  All  the  samples  that  have  survived  the  deposition 
were  found  in  the  predicted  safe  region  of  the  derived  design  space 
shown  in  Fig.  12.  More  speaking,  the  safe  fabricated  buckled  films  of 
the  ratios  of  side-length-to-thickness  (a/h)  were  exposed  to  the 
effective  compressive  stresses  values  ao  obtained  from  the  curva¬ 
ture  measurements  where  the  etching  effects  were  accounted  for 
through  Eq.  (14),  (see  Section  4.2,  Appendix  B  and  C)  and  the  cor¬ 
responding  points  (a/h,  ao),  marked  as  O,  were  all  found  above  the 
curves  Ci  and  C2  limiting  the  safety  regions  in  Fig.  12.  Moreover,  the 
8YSZ  samples  deposited  at  400  °C  under  2  mTorr  pressure  were  all 
cracked  and  actually  found  in  the  predicted  unsafe  region  where 
the  corresponding  points  (a/h,  ao)  are  marked  as  <g>. 

6.  Conclusion 

The  buckling  of  YSZ  thin-film  membranes  in  ^SOFC  is  analyzed 
by  applying  the  Rayleigh-Ritz  approach.  A  posteriori  analysis  of  the 
stress  measurements  is  suggested  to  derive  the  effective  residual 
stress  in  the  free-standing  membrane.  The  simulation  results  are 
found  to  be  in  agreement  with  the  experimental  measurements. 
The  thin  membrane  electrolyte  is  shown  to  survive  in  an  advanced 
stage  of  post-buckling  and  an  extended  area  of  safety  region  is 
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Fig.  12.  A  design  space  for  thin  film  deposition.  Safe  regions  are  denoted  as  Sj  where  unsafe  regions  (regions  of  to  avoid  from  design)  are  denoted  as  Yj,  i  =  1,2, 3, 4.  The  curves  Ci  and 
c2  are  the  thresholds  of  controlled  stresses  in  the  primary  and  the  secondary  post-buckling  regions.  These  curves  are  the  delimiters  of  the  safe  zones  Si  and  S2  of  controlled  stresses 
where  max  <  at  and  a0  >  <rc.  The  simulations  are  performed  with  Ritz  expansion  of  degrees  m  =  4  and  p  =  5,  for  Young’s  modulus  E  =  240  GPa,  Poisson  ratio  v  =  0.2  and 
tensile  strength  at  =  250  MPa.  All  the  8YSZ  samples  (square  films  of  300  nm  thickness  with  different  side-length  values)  deposited  at  different  temperatures  under  20  mTorr 
pressure  have  survived  the  deposition  and  were  found  in  the  predicted  safe  region  marked  as  O .  The  8YSZ  samples  deposited  at  400  °C  under  2  mTorr  pressure  were  all  cracked  and 
found  in  the  predicted  unsafe  region  marked  as  ® . 


shown  in  a  post-buckling  design  space.  An  advanced  exploitation  of 
the  energy  approach  would  be  the  prediction  of  post-buckling 
design  space  for  different  plate  geometries  (square,  rectangular, 
circular)  for  an  optimal  design  for  YSZ  manufacturing.  This  would 
be  the  topic  of  a  future  work.  The  method  used  in  this  work  is  us¬ 
able  as  a  step  in  the  study  of  the  effects  of  the  integration  of  elec¬ 
trodes  (that  have  a  different  thermal  stress  behavior)  with  the  YSZ 
electrolyte.  This  would  also  be  one  of  the  topics  of  a  new  paper  on 
the  prediction  of  the  safety  of  the  multi-layers. 
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Appendix  A.  Energy  storage  modes 

When  the  film  is  in  an  unbuckled  shape  for 
o0  >  =  -0.538  MPa,  the  potential  energy  is  stored  exclusively 

as  membrane  energy  which  is  “strongly  increasing”  as  the  in-plane 
residual  compression  is  augmented  in  pre-buckling  stage,  see 
Fig.  A.13(a).  A  film  in  primary  buckling  stage  exhibits  a  “lower  in¬ 
crease”  of  the  potential  energy.  Indeed  in  Fig.  A.13(a),  it  can  be 
observed  that  the  post-buckling  potential  energy  curve  has  the 
slope  of  the  tangent  to  the  pre-buckling  curve  at  the  transition 
point  of  first  buckling  i.e.  a0  =  oQ^  =  -0.538  MPa.  Flence,  buck¬ 
ling  allows  a  more  efficient  minimization  of  the  potential  energy  by 


an  exchange  from  membrane  energy  to  bending  energy.  The 
membrane  shear-energy  storage  is  minor  in  the  first  buckling  and 
this  corresponds  to  the  axis-symmetry  shape  of  the  buckled  film.  In 
secondary  buckling,  a  different  mode  of  energy  storage  is  taking 
place  with  a  considerable  amount  of  membrane  shear  energy 
Fig.  A.13(b),  this  can  only  be  ensured  by  breaking  the  mirror  sym¬ 
metry  of  the  deformation. 


-o-o  [MPa| 


(a) 


-0-o  [MPa] 


(b) 


Figure  A.13.  Distribution  of  potential  energy  in  different  stages  of  loading  a0 :  (a)  pre 
and  first  buckling  and  (b)  second  buckling.  Ritz  degree  m  =  4,  p  =  5. 
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Appendix  B.  Film-substrate  misfit  strain 
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Figure  B.14.  A  representation  of  stress  distribution  in  the  cross-section  of  the  sample  consisting  of  YSZ  film  deposited  on  substrate-supported  silicon  nitride  film  (a  sample  without 
etching). 


In  this  appendix,  the  misfit  strain  at  the  interface  between  film 
and  substrate  is  extracted  from  the  curvature  measurements  in  the 
case  of  the  substrate-supported  film.  This  estimation  will  be  used  in 
Appendix  C  as  interface  condition  to  obtain  the  effective  residual 
stress  in  the  case  of  free-standing  membrane.  The  stress  values  in  a 
case  of  substrate-supported  film  (without  etching)  is  obtained  by 
superposing  the  following  stress  contributions: 


_  al+abs(hs/2)+alNS 

tmis  fit  —  Bi  B2 


+ 


4hA 

hsB2) 


a 


C  _  °SNS 

f  B2  5 


(B.4) 


where  B\  denotes  the  biaxial  modulus  of  YSZ  supported  film. 


•  The  residual  stresses  in  Si3N4-Si-Si3N4  substrate,  Fig.  7(b),  which 
represent  a  misfit  interaction  between  Si3N4  and  Si  material. 
From  the  curvature  measurements  (see  Fig.  7,  steps  (a)  and  (b)) 
we  have  (j^NS  =  128  MPa  the  tensile  stress  in  the  Si3N4  film. 
Since  the  system  is  not  subjected  to  an  external  traction  the 
force  equilibrium  yields  a  compressive  stress  in  silicon  wafer: 

°SNS  ~  h3  °SNS- 

•  The  stresses  obtained  from  the  curvature  measurements  (see 
Fig.  7,  steps  (b)  and  (c))  which  are  compressive  in  YSZ  film  where 
af  =  -1350  MPa  and  tensile  in  Si3N4-Si-Si3N4  substrate  (deno¬ 
ted  by  subscript  “s”).  Based  on  the  force  equilibrium  we  obtain 


<4  =  (B-1) 

where  h\  is  the  YSZ  film  thickness  and  hs  =  2/12  +  h3  is  the  substrate 
thickness  with  h2  and  h3  are  the  Si3N4  film  coating  thickness  and 
the  wafer  Si  thickness  respectively. 

•  The  bending  stress  denoted  by  o-j?  in  Si3N4-Si-Si3N4,  it  is  written  as 


Appendix  C.  Residual  compressive  stress  in  the  free-standing 
YSZ  film 
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Figure  C.15.  A  representation  of  stress  in  the  free-standing  YSZ  and  Si3N4  membranes. 


Considering  the  free-standing  YSZ-Si3N4  double  layer,  see  Fig.  8, 
in  an  unbuckled  state  (this  corresponds,  however,  to  either  a  stable 
pre-buckling  state  under  low  in-plane  compressions,  r0  >  *  or  to 

an  assumed  unstable  equilibrium  state  of  unbuckled  film  under  high 
in-plane  compressions,  r0  <  see  paths  FP1  and  FP2  in  Fig.  5). 

The  interaction  between  the  two  layers  is  effected  by  the  forces 
and  the  moment  acting  at  the  interface.  Similar  to  the  approach 
introduced  by  Refs.  [17],  and  used  later  by  Refs.  [18],  [19]  and  [20], 
in  this  work  each  layer  is  modeled  as  a  beam  and  the  aforemen¬ 
tioned  interface  forces  and  moment  are  converted  to  equivalent 
stress  (or  force)  and  moment  system  referred  to  the  centroidal  axis 
of  each  layer.  The  bending  stresses  are: 


<*(*)  =  -j^z,  Vze 


hs  hs 

T’  ~2 


(B.2) 


where  R  is  the  measured  curvature  radius  and  Bs  is  the  average 
biaxial  modulus  of  the  substrate  Bs  =  2B2/12  +  B3h3lhs  whereas 
Bi  =  Bi/1  -  v\,  i  =  2,  3,  are  the  biaxial  modulus  of  Si3N4  film  and  Si 
wafer  respectively.  Further,  Stoney’s  relation  of  film  stress 

U2  D 

°f  =  ~6Rih  substituted  into  Eq.  (B.2)  giving: 


°YSz(.Z)  —  cT“Z’ 
Kfs 


hi  h T 
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(Cl) 
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(B-3) 


The  bending  stress  in  the  film  is  disregarded  because  of  its  low 
bending  stiffness  compared  to  the  thick  substrate.  Now  the  misfit 
strain  rmiS  fit  at  the  interface  between  YSZ  and  Si3N4  films  where 
z  =  hs! 2  is  obtain  as: 


where  a  plate  curvature  of  radius  RfS  (free-standing  curvature 
radius)  is  assumed. 

The  residual  compressive  stress  in  YSZ  membrane  (denoted  here 
as  0ySZ)  cannot  be  estimated  from  the  classical  Stoney’s  formula 
since  the  condition  of  thick  substrate  (h2  >>>  hi)  is  not  satisfied. 
Considering  the  torque  freedom  condition  at  the  middle  plane  of 
S^N4  layer  (Fig.  C.15),  the  moment  equilibrium  is  written  as: 
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Substituting  for  o\sz  and  <7!?^  from  Eqs  (C.1,  C.2)  into  Eq.  (C.3),  a 
new  formulation  (generalization)  of  Stoney’s  relation  is  obtained  as 
(Stoney-Safa): 


rC  ftfo+ftlBa 

YSZ  6Rfshi  (h]  +  h2) 


(C.4) 


For  a  very  thin  film  (hi  <<<  h2),  this  equation  gives  the  classical 
Stoney’s  formula.  The  tensile  stress  in  the  Si3N4  membrane  is  ob¬ 
tained  from  force  equilibrium  conditions: 


jS13N4 


(C.5) 


The  misfit  strain  £miSfit  at  the  interface  between  YSZ  and  Si3N4  is 
obtained  as  a  function  of  the  free-standing  membrane  (post¬ 
etching)  curvature  radius  RfS: 


^misfit 


^YSZ+^YSZ^) 


“BT 


A*) 


1  (  h3B^+h3B2  hi  h3Br+h3B2  h2\ 

~  RfS  y  6h1(h1+h2)B1  2  6h2(h,+h2)B2  2  )' 


(C.6) 


This  corresponds  to  the  relation  derived  elsewhere  see  e.g.  Refs. 
[21]  and  [22].  The  misfit  strain  between  the  two  membranes  is 
etching-independent,  and  then,  has  the  same  value  for  both  the 
substrate-supported  film  and  free  -standing  membrane  cases. 
Substituting  for  £misfit  from  Eq.  (C.6)  into  Eq.  (B.4),  we  obtain  an 
expression  for  the  curvature  radius  RfS  of  the  free-standing  mem¬ 
brane  for  an  assumed  unbuckled  state. 

Therefore,  the  compressive  stress  in  the  free-standing  YSZ  film 
is  obtained  by  substituting  the  aforementioned  obtained 


expression  for  RfS  in  Eq.  (C.4).  It  is  an  explicit  function  of  the 
measured  stresses  of  and  ofNS: 


(^_4fe)<7f  ~asns 
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